Quantum entanglement between the electron clouds of nucleic acids in DNA 
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We model the electron clouds of nucleic acids in DNA as a chain of coupled quantum harmonic 
oscillators with dipole-dipole interaction between nearest neighbours resulting in a van der Waals 
type bonding. Crucial parameters in our model are the distances between the acids and the cou- 
pling between them, which we estimate from numerical simulations [J. We show that for realistic 
parameters nearest neighbour entanglement is present even at room temperature. We quantify the 
amount of entanglement in terms of negativity and single base von Neumann entropy. We find 
that the strength of the single base von Neumann entropy depends on the neighbouring sites, thus 
questioning the notion of treating single bases as logically independent units. We derive an analyt- 
ical expression for the binding energy of the coupled chain in terms of entanglement and show the 
connection between entanglement and correlation energy, a quantity commonly used in quantum 
chemistry. 
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INTRODUCTION 



The precise value of energy levels is of crucial impor- 
tance for any kind of interaction in physics. This is also 
true for processes in biological systems. It has recently 
been shown for the photosynthesis complex FMO [2]-|5[ 
that maximum transport efficiency can only be achieved 
when the environment broadens the systems energy lev- 
els. Also for the olfactory sense the energy spectra of 
key molecules seem to have a more significant contribu- 
tion than their shape @. In Q the possibility of in- 
tramolecular refrigeration is discussed. A common theme 
of these works is the system's ability to use non-trivial 
quantum effects to optimise its energy levels. This leads 
to the question whether a molecule's energy levels are 
only determined by its own structure, or if the environ- 
ment shapes the molecule's energy level? Entanglement 
between system and environment is a necessary condi- 
tion to alter the system's state. Here we study the influ- 
ence of weak chemical bonds, such as intramolecular van 
der Waals interactions, on the energy level structure of 
DNA and discuss its connection to entanglement. To de- 
scribe the van der Waals forces between the nucleic acids 
in a single strand of DNA, we consider a chain of cou- 
pled quantum harmonic oscillators. Much work has been 
done investigating classical harmonic oscillators. How- 
ever, this cannot explain quantum features of non-local 
interactions. Also, classical systems can absorb energy 
quanta at any frequency, whereas quantum systems are 
restricted to absorb energy quanta matching their own 



energy levels. This is of importance for site specific DNA- 
Protein interaction, as the probability of a protein to bind 
to a specific sequence of sites in DNA is governed by the 
relative binding energy Q. 

Our work was motivated by a numerical study on the 
importance of dispersion energies in DNA [l| . Dispersion 
energies describe attractive van der Waals forces between 
non-permanent dipolcs. Recently their importance to 
stabilise macromolecules was realised [t| [Io[ . Modelling 
macromolecules, such as DNA, is a tedious and complex 
task. It is currently nearly impossible to fully quantum 
mechanically simulate the DNA. Quantum chemistry has 
developed several techniques that allow the simulation 
of DNA with simplified dynamics. In [l[ the authors 
first quantum mechanically optimise a small fragment of 
DNA in the water environment. Secondly, they "per- 
formed various molecular dynamics (MD) simulations in 
explicit water based either fully on the empirical poten- 
tial or on more accurate QM/ MM MD simulations. The 
molecular dynamics simulations were performed with an 
AMBER parm9916 empirical force field and the follow- 
ing modifications were introduced in the non-bonded part, 
which describes the potential energy of the system (see eq 
1 ) and is divided into the electrostatic and Lennard- Jones 
terms. The former term is modelled by the Coulomb in- 
teraction of atomic point- charges, whereas the latter de- 
scribes repulsion and dispersion energies, " 
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where the strength of the dispersion energy is scaled 
with the parameter e. For e = 1 the dynamics of 
the DNA strand is normal. For a weaker dispersion, 
e = 0.01, there is in increase of 27% in energy in the 
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DNA. This increase of energy induces the unravelling of 
the double helix to a flat, ladder-like DNA. Many factors 
contribute to the spatial geometry of DNA, e.g. water 
interaction, the phosphate backbone, etc. However, 
one of the strongest contributions is the energy of the 
electronic degree of freedom inside a DNA strand, which 
is well shielded from interactions with water. Stronger 
interaction (e = 1) allows the electrons clouds to achieve 
spatial configurations that require less structural energy. 
This allows a denser packing of the electron charges 
inside the double helix. 

Here we investigate with a simple model of DNA 
whether continuos variable entanglement can be present 
at room temperature, and how this entanglement is 
connected to the energy of the molecule. There are many 
technically advanced quantum chemically calculations 
for van der Waals type interaction, i.e. 11] . The aim of 
this work is not to provide an accurate model, but to 
understand underlying quantum mechanical features and 
their role in this biological system. Also, there are many 
parallel developments between quantum information 
and quantum chemistry. This work bridges the concepts 
of entanglement and dispersion energies between the two 
fields. Finally, the advantages of quantifying chemical 
bonds in terms of entanglement were already mentioned 
in [ljj ■ Here we give the first example of a system whose 
chemical bonds are described by entanglement. 



DISPERSION ENERGIES BETWEEN NUCLEIC 
ACIDS 

The nucleic bases adenine, guanine, cytosine and 
thymine are planar molecules surrounded by ir electron 
clouds. We model each base as an immobile positively 
charged centre while the electron cloud is free to move 
around its equilibrium position, see Fig. [TJ There is no 
permanent dipole moment, while any displacement of the 
electron cloud creates a non-permanent dipole moment. 
Denoting the displacement of two centres by (x, y, z), we 
assume the deviation out of equilibrium |(x, y, z)| to be 
small compared to the distance r between neighboring 
bases in chain. The displacement of each electron cloud 
is approximated to second order and described by a har- 
monic oscillator with trapping potential fl that quanti- 
fies the Coulomb attraction of the cloud to the positively 
charged centre. A single DNA strand resembles a chain 
of harmonic oscillators, see Fig. [2j where each two neigh- 
boring bases with distance r have dipolc-dipolc interac- 
tion. 

The Hamiltonian for the DNA strand of N bases if 
given by 





FIG. 1: This graphic shows a sketch of a DNA nucleic acid. 
The mostly planar molecules are divided into the positively 
charged molecule core (red) and the negatively charged outer 
7r electron cloud (blue-yellow). In equilibrium the centre of 
both parts coincide, thus there is no permanent dipole. If the 
electron cloud oscillates around the core, a non permanent 
dipole is created [TJ- The deviation out of equilibrium is de- 
noted by (ar, y, z). The corresponding dipole is jl — Q(x, y, z). 
This oscillation might be caused by an external field, or in- 
duced by quantum fluctuations, as it is given in a DNA strand. 



where d denotes the dimensional degree of freedom, and 
the dipole potential 

r- 1 
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with flj — Q(xj, yj, Zj) dipole vector of of site j and fjv 
normalised distance vector between site j and j + Due 
to symmetry fjv is independent of j. The dimensionless 
scaling factor e is varied to study the effects on entan- 
glement and energy identical as in [l|. In order to com- 
pare our model with [H , we consider 'normal' interaction, 
where the dipole-dipole interaction has full strength mod- 
elled by e = 1 and 'scaled' interaction, where the dipole- 
dipole interaction is reduced to a hundredth of the origi- 
nal strength modelled by e = 0.01. The distance between 
neighbouring bases in DNA is approximately ro = 4.5A. 
For generality we will not fix the distance. 
In general the single strand of DNA will not be per- 
fectly linear and thus the dipole potential has coupling 
terms of the form xz etc. Detailed analysis following [14| 
shows that the energy contribution from the cross cou- 
pling terms is small, and they will be ignored here. This 
leads to the interaction term 
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The different signs for x,y and z reflect the orientation 
of the chain along z direction. 

A discrete Fourier transformation of the form 
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FIG. 2: This graphic shows a sketch of a single DNA strand. 
The chain is along z direction. Each bar in the single strand 
DNA represents one nucleic acid: adenine, thymine, guanine 
or cytosine. Around the core of atoms is the blue outer elec- 
tron cloud. The oscillation of these electron clouds is modelled 
here as non-permanent harmonic dipoles, depicted by the ar- 
rows, with trapping potential fid in dimension d — x,y, z. 



decouples the system into independent phonon modes. 
These modes can be diagonalized by introducing creation 
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a d , . This results in the dispersion relations 
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where n^.i — a\ jftd.z is the number operator of mode / in 
direction d. 



The trapping potentials fid can be linked to experimen- 
tal data (see table HJ) through the relation fid = \j 
, where ay is the polarizability of the nucleid base. So far 
we did not discuss the number of electrons in the cloud. 
Both the trapping potential fl 2 d as well as the interaction 

term ^— depend linearly on the number of electrons, and 
thus the dispersion frequencies w 2 , ; have the same de- 
pendance. The quantities of interest of this paper are 
entanglement and energy ratios, which are both given by 
ratios of different dispersion frequencies and are thus in- 
variant of the number of electrons involved. In Table UU 
we assumed the number of interacting electrons to be 
one, but our final results arc independent of this special 
choice. 

TABLE I: Numerical values for polarizability of different nu- 
cleid acid bases [TJ|] in units of lau = 0.164 ■ lCT 40 Fm 2 . 
The trapping frequencies are calculated using the formula 

n = 



— — and are given in units of 10 is Hz. 



nucleic acid 


Qi 


a y 






fly 


£l z 


adenine 


102.5 


114.0 


49.6 


4.1 


3.9 


6.0 


cytosine 


78.8 


107.1 


44.2 


4.7 


4.1 


6.3 


guanine 


108.7 


124.8 


51.2 


4.0 


3.8 


5.9 


thymine 


80.7 


101.7 


45.9 


4.7 


4.2 


6.2 



Although the values for the four bases differ, all show 
similar fl x w tt y (transverse), while there is an increase 
of 50% in the longitudinal direction, Cl z sa ^^x,y In 
the following we will approximate the chain to have the 
same value of trapping potential at each base. In x,y 
direction we will use fl x>y = 4- 10 15 Hz, and in z direction 
Q z = 6 • l0 15 Hz. 



ENTANGLEMENT AND ENERGY 

We now clarify the influence of entanglement on 
energy. We will also derive an analytic expressions for 
the change in binding energy depending on entanglement 
witnesses. 

The chain of coupled harmonic oscillator is entan- 
gled at zero temperature, but is it possible to have 
entanglement at room temperature? There is a 
convenient way to calculate a criterion for nearest 
neighbour entanglement for harmonic chains [l6j |. 
which compares the temperature T with the coupling 
strength w between neighbouring sites. In general, for 
2 ^ T < 1 one can expect entanglement to exist. Here 
the coupling between neighbouring clouds is given by 

i e 1 / 4 !^ • 10 15 Hz for r = 4.5l, which 

0.05 for e = 1 and 0.16 for e = 0.01. 
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This means that the coupling between electron clouds 
is dominant compared to the temperature, and thus 
implies the existence of entanglement even at biological 
temperatures. An exact method to quantify the amount 
of entanglement in harmonic states it the violation of 
one of the two inequalities, related to the covariance 
matrix the state 
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< Si = ¥ ((dj + d ]+l f) ((p dd - Pw f) - 1 (8) 
< S 2 = ^ ((dj - d j+l f) (ip d!j + PW ) 2 ) - 1 (9) 



with gL- 



j position operator of site j in direction d and 
Pd.j corresponding momentum operator. If one of the 
inequalities is violated, the sites j and j + 1 are en- 
tangled. The negativity, a widely used measure for 
entanglement, is calculated using the formula Neg = 
5Zfc=i max [0, — In \/Sk + l] • The negativity measures 
the amount of entanglement between two subsystems. It 
can be directly calculated from space and momentum op- 
erator expectation values, namely the above defined Si ,2 
criteria. The amount of negativity between neighbouring 
bases for room temperature is shown in Fig. [31 For the 
normal coupling there is substantially more entanglement 
present than for the scaled interaction. This correlates 
with the amount of binding energy found in [![, where 
the DNA with normal coupling has a lower energy than 
the DNA with scaled coupling. 

The above result motivates the question whether the 
binding energy can be expressed in terms of entanglement 
measures. In the limit of long distances, an analytical 
expression connects the amount of binding energy in the 
chain of oscillators with the values of S\ : 2- Due to the 
strong coupling the chain of oscillators is effectively in 
its ground state, which we will assume in the following 
analysis. 

The dispersion relations of the electron cloud oscilla- 
tions can be expanded for large distances, i.e. r 
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Inserting this expansion into the entanglement criterion 
5-2 gives: 
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while the corresponding expression for S^i has a positive 
value. A similar expansion of the dispersion relation in 
x direction leads to: 
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FIG. 3: This graphic shows the nearest neighbour negativity 
as a function of distance between sites in A at T = 300A'. 
The three upper curves are for scaling factor e = 1, the lower 
two curves are for scaling factor e = 0.01. The red curve is for 
z direction and fl z = 6 ■ 10 15 Hz. The blue and green curve 
are for x direction and Q x = 4 • W 15 Hz and Q, x = 3 • W 15 Hz. 
The negativity for e = 0.01 is much smaller than in the un- 
sealed case. The amount of negativity strongly depends on 
the distance r between sites and the value of trapping poten- 
tial Q. The lower the potential, the higher the negativity. A 
typical distance between neighbouring base pairs in DNA is 
approximately r — 4.5A. Along the chain (^-direction) the Si 
criterion is violated, whereas transversal to the chain S2 (x- 
direction) is violated. This reflects the geometry of the chain. 
Along the main axes of the chain energy is reduced by corre- 
lated movement. Transversal to the chain it is energetically 
better to be anti-correlated. 



This implies that nearest neighbor (n.n.) electronic 
clouds are entangled even at large distances. However 
the amount of entanglement decays very fast. We will 
now compare this result with the binding energy in the 
ground state. The binding energy is defined as the dif- 
ference of energy of the entangled ground state and any 
hypothetical separable configuration 
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(14) 

This definition is analogous to the definition of correla- 
tion energy in chemistry [l8T ] . The first approximation to 
the full Schrodinger equation is the Hartree-Fock equa- 
tion and assumes that each electron moves independent 
of the others. Each of the electrons feels the presence of 
an average field made up by the other electrons. Then 
the electron orbitals are antisymmetrised. This mean 
field approach gives rise to a separable state, as antisym- 
metrisation does not create entanglement. The Hartrcc- 
Fock energy is larger than the energy of the exact solution 
of the Schrodinger equation. The difference between the 
exact energy and the Hartree-Fock energy is called the 



5 



correlation energy 

Ecorr = E exact — EhF ■ (15) 

Our definition of binding energy is a special case of 
the correlation energy, but we restrict our analysis here 
to phonons (bosons) instead of electrons. Our model 
describes the motional degree of freedom of electrons, 
namely the displacement of electron clouds out of equi- 
librium. We show for this special case that the amount of 
correlation energy is identical to entanglement measures. 
Expanding the binding energy for r 3 — > oo, the leading 
term is of order 




(16) 

since the first order vanishes due to symmetry and simi- 
larly for x direction: 

E x ,bind ~ 7T—Si . (17) 

o 

Eq. [THl [T7] show a simple relation between the entangle- 
ment witnesses Sx 2 and the binding energy of the chain 
of coupled harmonic oscillators. The stronger the en- 
tanglement, the more binding energy the molecule has. 
Interestingly, along the chain the S\ criterion is violated, 
whereas transversal to the chain S2 is violated. This re- 
flects the geometry of the chain. Along the main axes 
of the chain energy is reduced by correlated movement. 
Transversal to the chain it is energetically better to be 
anti-correlated. This means that the entanglement wit- 
nesses Si t 2 not only measure the amount of binding en- 
ergy, but also the nature of correlation which gives rise to 
the energy reduction. This relation motivates the search 
for entanglement measures describing the binding ener- 
gies of complex molecules. While the binding energy just 
measures energy differences the corresponding entangle- 
ment measures reflect more information. Without corre- 
lations between subsystems there would not be a chem- 
ical bond. It is precisely the purpose of entanglement 
measures not only to quantify, but also to characterise 
these correlations. 

APERIODIC POTENTIALS AND 
INFORMATION PROCESSING IN DNA 

In the above calculations we assumed a periodic po- 
tential, which allowed us to derive analytical solutions. 
Here we investigate the influence of aperiodic potentials 
and discuss the robustness of the previous conclusions. 
Firstly we note that the potentials for different nucleic 
acids do not differ significantly, see table HI Hence one 
would intuitively assume that a sequence of different lo- 
cal potentials changes the amount of entanglement but 



does not destroy it. To check this intuition more thor- 
oughly, one can use the phonon frequencies of the ape- 
riodic chain of oscillators. For a finite one-dimensional 
chain of 50 bases without periodic boundary conditions 
and with the sequence of nucleic acids randomly chosen, 
we solve the resulting coupling matrix numerically. The 
smallest dispersion frequency determines the thermal ro- 
bustness; the smaller the frequencies w; the larger the 
probability that the thermal heat bath can excite the 
system. Sampling over 1000 randomly chosen sequences 
yielded min(a;;) = 3.2 ■ 10 15 Hz as smallest dispersion fre- 
quency. Comparing this with the thermal energy gives 

2kB n°° K ~ °' 03 ' wnich is sti11 ver y smalL 

Thus the thermal energy is more than 20 times smaller 

than the smallest phonon frequency, which allows us to 
continue working with the ground state of the system. 
Different sequences will cause fluctuations in the amount 
of entanglement in the chain of bases. We determine for 
each string the average of single site von Neumann en- 
tropy and compare it with the classical amount of infor- 
mation measured by the Shannon entropy of each string. 
The von Neumann entropy of a single site j is obtained 
following [l?! with the formula 

MrJ-2±iM(2±i)-2fi,.(2^) (18, 

where Tj = j-^J ' (x^)(p Xj ) : is the symplcctic eigenvalue of 

the covariance matnrix of the reduced state. 
To check whether the relative frequency of A,C,G 
and T influences the amount of entanglement within 
the coupled chain of oscillators, we also calculate the 
classical Shannon entropy of each string. Fig. @] shows 
the average amount of single site quantum entropy vs. 
classical entropy. There is, within this model, no direct 
correlation between classical and quantum entropy. For 
the same amount of Shannon entropy, i.e. same relative 
frequencies of A,C,G and T, the value of quantum cor- 
relations varies strongly between around vNE = 0.007 
and vNE = 0.025. We note that for achieving a 
comparable amount of local disorder by thermal mixing 
a temperature of more than 2000i^ is needed. This is a 
quantum effect without classical counterpart. Each base 
without coupling to neighbours would be in its ground 
state, as thermal energy is small compared to the energy 
spacing of the oscillators. As the coupling increases, 
the chain of bases evolves from a separable ground 
state to an entangled ground state. As a consequence 
of the global entanglement, each base becomes locally 
mixed. This feature cannot be reproduced by a classical 
description of vibrations. When a classical system is 
globally in the ground state, also each individual unit is 
in its ground state. Although it is already well known 
that globally entangled states are locally mixed, little 
is known about possible consequences for biological 
systems. In the following paragraph we discuss one such 
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FIG. 4: This graphic shows the average single site von Neu- 
mann entropy of a chain of nucleic acids dependant on the 
classical Shannon entropy of the string. Each string contains 
50 bases with a random sequence of A,C,G, or T. The dis- 
tribution of nucleic acids determines the classical Shannon 
entropy. For each nucleic acid we used the value of polariz- 
ability of Table [I] in x direction. The distance between sites is 
r — 4.5 • lCP 10 m. The plot has a sample size of 1000 strings. 
There is no direct correlation between quantum and classical 
information. The average amount of von Neumann entropy 
varies strongly for different sequences. 



we consider it more accurate to describe the processing 
of genetic information by quantum channels, as the in- 
teractions between molecules are determined by laws of 
quantum mechanics. 



TABLE II: Numerical values for the von Neumann entropy 
of site 9 (Adenine) in a chain with open boundary condition 
containing 17 bases. The bases 1-7 and 11-17 are taken to be 
Adenine. The column gives the nucleic acid of site 8, the rows 
of site 10. The von Neumann entropy of site 9 varies with its 
neighbours. 





Adenine 


Cytosine 


Guanine 


Thymine 


Adenine 


0.077 


0.082 


0.078 


0.081 


Cytosine 


0.082 


0.079 


0.084 


0.078 


Guanine 


0.078 


0.084 


0.079 


0.083 


Thymine 


0.081 


0.078 


0.083 


0.078 



CONCLUSIONS AND DISCUSSION 



quantum effect on the information flow in DNA. 



How much information about the neighbouring sites is 
contained in the quantum degree of freedom of a single 
base? Is it accurate to describe a single nucleic acid as an 
individual unit or do the quantum correlations between 
bases require a combined approach of sequences of nucleic 
acids? The single site von Neumann entropy measures 
how strongly a single site is entangled with the rest of 
the chain and is therefore a good measure to answer this 
question. In the following we considered a string with 
17 sites of a single strand DNA. Site 9 as well as sites 
1-7 and 11-17 are taken to be Adenine. The identity of 
nucleic acids at sites 8 and 10 varies. Table lU shows 
the resulting von Neumann entropy of site 9 dependent 
on its neighbours. The value of a single site depends on 
the direct neighbourhood. There is, for example, a dis- 
tinct difference if an Adenine is surrounded by Cytosine 
and Thymine (vNE = 0.078) or by Cytosine and Gua- 
nine (vNE = 0.084). On the other hand, in this model 
there is little difference between Adenine and Guanine in 
site 8 and Guanine in site 10. Of course this model has 
not enough precision to realistically quantify how much 
a single site knows about its surroundings. Nevertheless 
it indicates that a single base should not be treated as an 
individual unit. When quantifying the information con- 
tent and error channels of genetic information, the anal- 
ysis is usually restricted to classical information trans- 
mitted through classical channels. While we agree that 
the genetic information is stored using classical informa- 
tion, e.g. represented by the set of molecules (A,C,G,T), 



In this paper we modelled the electron clouds of nu- 
cleic acids in a single strand of DNA as a chain of coupled 
quantum harmonic oscillators with dipolc-dipole interac- 
tion between nearest neighbours. Our main result is that 
the entanglement contained in the chain coincides with 
the binding energy of the molecule. We derived in the 
limit of long distances and periodic potentials analytic 
expressions linking the entanglement witnesses to the en- 
ergy reduction due to the quantum entanglement in the 
electron clouds. Motivated by this result we propose to 
use entanglement measures to quantify correlation en- 
ergy, a quantity commonly used in quantum chemistry. 
As the interaction energy given by hio is roughly 20 times 
larger than the thermal energy IcbSOOK the motional 
electronic degree of freedom is effectively in the ground 
state. Thus the entanglement persists even at room tem- 
perature. Additionally, we investigated the entanglement 
properties of aperiodic potentials. For randomly chosen 
sequences of A,C,G, or T we calculated the average von 
Neumann entropy. There exists no direct correlation be- 
tween the classical information of the sequence and its 
average quantum information. The average amount of 
von Neumann entropy varies strongly, even among se- 
quences having the same Shannon entropy. Finally we 
showed that a single base contains information about its 
neighbour, questioning the notion of treating individual 
DNA bases as independent bits of information. 
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